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Abstract 

We compare the physics potential of the upcoming neutrino oscillation experiments 
Daya Bay, Double Chooz, NO^A, RENO, and T2K based on their anticipated nominal lu- 
minosities and schedules. After discussing the sensitivity to #13 and the leading atmospheric 
parameters, we demonstrate that leptonic CP violation will hardly be measurable without 
upgrades of the T2K and NOz/A proton drivers, even if #13 is large. In the presence of the 
proton drivers, the fast track to hints for CP violation requires communication between 
the T2K and NOz/A collaborations in terms of a mutual synchronization of their neutrino- 
antineutrino run plans. Even in that case, upgrades will only discover CP violation in a 
relatively small part of the parameter space at the 3a confidence level, while 90% confidence 
level hints will most likely be obtained. Therefore, we conclude that a new facility will be 
required if the goal is to obtain a significant result with high probability. 
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1 Introduction 



Neutrino oscillations have been firmly established in the last ten years or so by a beautiful 
series of experiments with neutrinos from the sun [1-6], the Earth's atmosphere [7,8], nu- 
clear reactors [9,10], and accelerators [11,12]. While these measurements have discovered 
and confirmed the dominant effective 2-flavor oscillation modes, it will be the purpose of 
the upcoming generation of experiments to discover sub-leading effects. This includes the 
following tasks: 

1. Determination of the small lepton mixing angle #13, 

2. Establishing CP violation (CPV) in neutrino oscillations for a value of the Dirac CP 
phase 5 cp 7^ 0, n, 

3. Identification of the type of the neutrino mass hierarchy (MH), which can be normal 
(Am^ > 0) or inverted (Aro^ < 0). 

There are several neutrino oscillation experiments currently under construction, which 
are expected to start data taking soon. These are the reactor neutrino experiments Dou- 
ble Chooz [13], Daya Bay [14], RENO [15] and the accelerator experiments T2K [16] and 
NOz/A [17]. The primary goal for all of these experiments is the discovery of the yet un- 
known mixing angle #13. In this work, we will revisit the nominal sensitivities of these 
projects to #13 taking care of the different nature of the experiments (z/ e disappearance in 
reactors versus z/ M — > z/ e or z/ M — > z/ e appearance in accelerators), and estimate the time evo- 
lution of the global discovery reach for 6> 13 based on the official schedules of the experiments. 
This analysis updates previous works [18-22] with respect to the now settled parameters for 
the considered experiments. Let us mention that eventually already the currently running 
experiments MINOS [12,23] and OPERA [24] might give a first hint for a non-zero value of 
#13 in case it is close to the present bound. We do not consider these experiments here, but 
focus on the above facilities currently under preparation, since they will clarify such hints 
with high significance. 

Furthermore, we will study if there is any chance to address items 2 and 3 above (CPV 
and MH) already with this set of upcoming experiments in case of a soon coming positive 
signal for 13 , which implies that #13 is relatively "large". It turns out that even in this 
most favorable case the sensitivity to CPV and MH of these experiments in their nominal 
configuration is marginal. Therefore we will explore the potential of "minor upgrades" to the 
proposed setups of T2K and NOz/A, based upon mostly existing equipment. These include 
a longer running time and an upgraded beam power for both experiments, and the addition 
of antineutrino running in T2K. Furthermore, we investigate the optimization potential of a 
coordinated combination of T2K, NOz/A, and Daya Bay. The purpose of this analysis is to 
estimate whether there will be any chance to have information on CPV or the MH around 
2022-2025, given the rather optimistic assumptions concerning beam performances and the 
size of #13. 

Such considerations are relevant in view of the current international effort towards a 
subsequent high precision neutrino facility [25,26]. These studies aim for a decision point 
concerning a future facility around 2012, where improved information on # 13 should be 
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Table 1: Summary of the standard setups at their nominal luminosities. 



available from the above mentioned experiments. If a non-zero value is established until 
then, an eminent question will be wether CPV and MH can be found by upgrades of the 
existing experiments, or indeed a new facility will be necessary. 

The outline of the paper is as follows. In Sec. El we describe the planned experiments 
and give some details on our simulations. In Sec. [3J we consider the nominal configurations 
of the experiments based on the information the experimental collaborations provide. We 
study the sensitivity to #i 3 (discovery potential as well as the case of large #i 3 ) and the im- 
provements to be expected in the leading atmospheric parameters #23 and [Aml-J. In Sec. IU 
we discuss the prospective time evolution of the sensitivity to 6^3 within the upcoming years. 
In Sec. [5j we consider possible upgrades for T2K and NOt/A focusing on the measurements 
of CPV and MH in the case of relatively large #13. Note that we arranged the Sees. [3] to [5] 
in an order from the most to the least established in terms of the data provided by the 
experimental collaborations. Summary and conclusions follow in Sec. [6l In the technical 
appendix we give details on our neutrino/antineutrino optimization algorithm. 

2 Experiment descriptions and simulation methods 

Below we describe in some detail the considered experimental setups and our simulation 
which we perform by using the GLoBES software [27,28]. The corresponding gib-files are 
available at the GLoBES web-page [27] including detailed technical information on the sim- 
ulation. In all cases our strategy is to follow as close as possible the original Letters of 
Intent (LOIs) or Technical Design Reports (TDRs). We have made sure that our sensitivi- 
ties agree with the "official" curves from the corresponding collaborations under the same 
assumptions. Table [1] summarizes the key parameters of the considered experiments. 

Reactor experiments look for the disappearance of v e , governed by #13, where the neu- 
trinos are produced in the nuclear fission processes in commercial nuclear power plants. All 
three experiments, Double Chooz, Daya Bay, and RENO, use a liquid scintillator doped 
with Gadolinium in order to exploit the coincidence of a positron and a neutron from the 
reaction u e + p — > e + + n. The two crucial parameters determining the final sensitivity are 
the total exposure (proportional to the detector mass times the time-integrated thermal 
power of the reactor over the lifetime of the experiment) and the systematic uncertainty, 
see, e.g., Ref. [20] for a discussion. All three proposals use the concept of near and far de- 
tectors in order to reduce uncertainties on the initial neutrino flux. Our implementation of 
systematic uncertainties is similar to the one of Ref. [29] . While for Double Chooz and Daya 
Bay detailed proposals are available, Refs. [13] and [14], respectively, the RENO experiment 
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is somewhat less documented and some properties have to be extrapolated from the other 
two proposals. The characteristics of the three reactor experiments are (see also Ref. [30] 
for a comparison study): 

Double Chooz: The Chooz power plant in France consists of two reactors with 4.3 GW 
thermal power each. There will be two detectors with 8.3 t fiducial mass each, a far 
detector at a distance of 1.0 km and 1.1 km from the two cores, and a near detector 
at a distance of 470 m and 350 m, respectively. Including efficiencies of 80%, detector 
dead times, and a load factor of 78% for the reactors, the event rates per year are 
8 x 10 4 and 1.5 x 10 4 for the near and far detectors, respectively. The uncertainty on 
the relative normalization of the detectors is assumed to be 0.6%. 

Daya Bay: The Chinese Daya Bay reactor complex consists currently of two pairs of reac- 
tor cores (Daya Bay and Ling Ao), separated by about 1.1 km. The complex generates 
11.6 GW of thermal power; this will increase to 17.4 GW by early 2011 when a third 
pair of reactor cores (Ling Ao II) is put into operation. For the 3 years nominal 
exposure we assume that all the mentioned six cores are operational. In total, eight 
detector modules will be installed, 20 t fiducial mass each. The far site consists of 4 
modules with distances from the three reactors of 1.985, 1.613, and 1.618 km. There 
will also be two near labs, each containing two detector modules at distances from 
the reactors between 360 and 530 m. The event rate per year is about 4 x 10 5 at the 
near sites and about 10 5 at the far site. The detector normalization is assumed to be 
0.18% uncorrelated between the detector modules. This corresponds to the final goal 
in terms of systematics without swaping detector modules. 

RENO: The RENO experiment will be located at the Yonggwang reactor complex in South 
Korea, which consists of 6 cores equally distributed along a straight line of about 
1.5 km, with a total thermal power of 16.4 GW. Two detectors with about 15 t 
fiducial mass each will be installed roughly at distances of 320 m and 1.4 km from the 
reactor line. Event rates per year will be about 6 x 10 5 at the near site and about 
3 x 10 4 at the far site. The relative detector normalization is assumed at 0.5%. 

The T2K and NOz/A superbeam experiments search for the appearance of electron neu- 
trinos in a beam of mainly muon neutrinos, from the decay of pions and kaons produced 
at a proton accelerator. Such beams unavoidably contain some intrinsic electron neutrinos 
which consist a background for the appearance search. Both experiments explore the off- 
axis technique to suppress the electron neutrino as well as the neutral current backgrounds 
and to achieve a more peaked beam spectrum. For both experiments, we also include the 
disappearance channels. 

T2K: The T2K experiment [16] sends a neutrino beam from the J-PARC accelerator to 
the Super-Kamiokande water Cerenkov detector with a fiducial mass of 22.5 kt at a 
distance of 295 km. Our simulation is based on publicly accessible sources as of 2008. 
We assume a 2.5 degree off-axis beam corresponding to 0.75 MW of beam power for 
5 years of neutrino running. The analysis of the v e appearance follows the thesis 
of M. Fechner [31]. Event rates and spectra for total signal and total background 
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match Fig. 6.8 of Ref. [31], and systematics are taken from chapter 6 of that reference. 
Our sensitivity from the appearance channel agrees with Fig. 6.17 of Ref. [31]. The 
disappearance analysis is taken from a talk by I. Kato given at Neutrino 2008 [32]. 
From that source we obtain rates and spectra for signal and non-QE background. 

NOfA: The description of NOz/A concerning the v e jv e appearance signals follows the 
proposal as of March 15, 2005 [17], the description of the disappearance signal is 
taken from Ref. [33]. We calibrated our event simulation to the numbers for signal, 
background, and efficiencies given in the October 2007 TDR [34] (Tabs. 6.2-6.4). 
NOi^A is a Totally Active Scintillator Detector (TASD) with a mass of 15 kt and is 
located at a distance of 810 km from NUMI beam source at Fermilab. The nominal 
luminosity is 18 x 10 20 protons on target for neutrinos and antineutrinos, each. This 
number of protons on target is assumed to correspond to 3 years running at 0.7 MW 
beam power. While this equal exposure to neutrinos and antineutrinos is the default 
assumption, we will also consider the sensitivities in case of neutrino running only 
in some of the next sections. We adopt systematical errors of 5% on the signal (not 
stated in the TDR) and 10% on the background (as given in the TDR). We use an 
analysis energy window from 1 to 3 GeV and assume the so-called medium energy 
(ME) beam configuration. 

For the sensitivity analyses we use the oscillation parameter values from Ref. [35]: 
Am 2 21 = 7.65 • 10" 5 eV 2 , |Am^| = 2.40 ■ 10" 5 eV 2 , sin 2 12 = 0.304, and sin 2 8 23 = 0.500, 
unless stated otherwise. We impose external la errors on Am^ (3%) and 9±2 (4%) as con- 
servative estimates for the current measurement errors [35], as well as Am 2 x (5%) for reactor 
experiments if analyzed without beam experiments. In addition, we include a 2% matter 
density uncertainty. 

3 Physics potential at nominal luminosities 

In this section, we discuss the physics potential of the experiments listed in Table [1] based 
on their nominal luminosities, i.e., the running times, detectors masses, and total numbers 
of protons (or thermal reactor powers) anticipated by the respective collaborations. 

3.1 Discovery potentials for 13 , MH, and CPV 

We show the discovery potentials for #13, MH, and CPV in Fig.[IJ These discovery potentials 
quantify for any given (true) sin 2 2# 13 for which fraction of possible (true) values of <5cp 
the corresponding quantity will be discovered^] For # 13 this means excluding the value 
#13 = 0. The mass hierarchy discovery potential is defined as the ability to exclude any 
degenerate solution for the wrong (fit) hierarchy at the chosen confidence level. Similarly, 
the CP violation discovery potential refers to the ability to exclude the CP conserving 
solutions 5qp = and Sep = vr for any degenerate solution at the chosen confidence level. 
These discovery potentials of course depend on the true values of 9 13 and 6qp and the true 

In this work we use the terms true value and simulated value synonymously in order to denote the values 
of parameters adopted to generate "data" , in contrast to the fit values used to fit these data. 
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Figure 1 : #13 , MH, and CPV discovery potential as fraction of true <5cp as a function of the true sin 2 2f?i3 
for the normal hierarchy (upper row) and inverted hierarchy (lower row) at the 90% CL. Note the different 
vertical scales in the different panels. 

hierarchy. In Fig. [T] we show for a given true value of sin 2 2# 13 (horizontal axis) and a given 
true hierarchy (upper row normal, lower row inverted) the fraction of all possible true values 
of (5cp for which the discovery can be achieved at the 90% confidence level. Hence, a fraction 
of (5cp of unity (or 100%) for a given sin 2 29i 3 corresponds to a discovery for any possible 
value of <5cp. 

The #i3 discovery potential (c/., left panels of Fig. [Q of the reactor experiments does 
not depend on <5cp since by convention this phase does not appear in the disappearance 
probability P ee . Furthermore, the probability is given to good approximation by an effective 
2-flavor expression: -P e r | act ~ 1 — sin 2 2t9 13 sin 2 (Am 2 1 L/4£'). Thanks to the large exposure, 
Daya Bay will have the best discovery potential among the reactor experiments of sin 2 2^13 = 
0.0066 at the 90% CL, compared to 0.018 for RENO and 0.033 for Double ChoozJJ In 
contrast, the — > v e appearance probability relevant for the beam experiments shows a 

2 Let us mention that the Daya Bay assumptions of a systematical error of 0.18%, fully uncorrelated 
among all detectors is more aggressive than for other reactor experiments. For example, if the systematic 
error is at the level of 0.6%, such as assumed in Double Chooz, the Daya Bay sensitivity of sin 2 2f5 13 = 0.0066 
deteriorates to sin 2 2f5 13 ~ 0.01. If on the other hand the systematic error is 0.38% and assumed to be fully 
correlated among modules at one site the limit would sin 2 2f5 13 ~ 0.012 [36]. See also the discussion in 
Ref. [30]. 
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dependence on the CP phase due to an interference term of the oscillations with Am^ and 
Am^. A discussion of the complementarity of reactor and beam experiments can be found, 
e.g., in Refs. [18,20,22]. Whether the best 6*i 3 discovery potential is obtained from Daya Bay 
or from one of the beams depends on the true Sep. For instance, for the normal hierarchy, 
T2K and NOz/A will have a discovery potential for slightly smaller values of sin 2 2#i3 than 
Daya Bay for about 30% of all possible 5qp values. 

The mass hierarchy measurement (cf., middle panels of Fig. [T]) requires NOz/A because 
of its relatively long baseline and therefore significant matter effects. T2K does not have 
any mass hierarchy discovery potential. For very large sin 2 2#i3, the mass hierarchy will be 
discovered for about 40-50% of all values of 5qp at 90% CL. There is no sensitivity left at 
this CL for sin 2 2#i3 < 0.04. Only minor improvement can be achieved if other experiments 
are added to NOz/A; the combination with T2K helps somewhat for large #13, whereas the 
reactor experiments contribute somewhat for small #13. 

For the discovery of CP violation (cf., right panels of Fig. [I]), neither T2K nor NOz/A 
alone have a substantial potential. Only the combination of these two experiments can 
measure CPV at 90% CL for up to 30% of all values of 5cp if the hierarchy is inverted. This 
hierarchy dependence appears because the matter effect in NOz/A leads to a more balanced 
neutrino-antineutrino statistics for the inverted hierarchy due to an matter enhancement of 
the antineutrino event rate. In the presence of the reactor experiments, however, the CP 
violation discovery potential improves significantly, especially for the normal hierarchy. If 
the hierarchy is normal, Daya Bay may in fact be the key prerequisite to obtain an early 
hint for CP violation. 

We stress that these sensitivities for MH and CPV are at the 90% CL, which by no means 
can be considered as a discovery. Increasing the CL leads to drastic loss in sensitivity and 
almost nothing can be said about CPV and MH at the 3a CL. 

3.2 The case of 13 just around the corner 

We show in Fig. [2] (normal simulated hierarchy) and Fig. [3] (inverted simulated hierarchy) 
how typical fits in the 6*i3-5cp plane would look like if 6*13 was large (sin 2 29 1 3 = 0.1) and Sep 
was close to maximal CP violation 5qp = vr/2 (upper rows) and <5cp = 37r/2 (lower rows). 
The fit contours for the right fit hierarchy are shaded (colored), the ones for the wrong fit 
hierarchy fit are shown as curvesJl 

The figures show the characteristics of the different classes of experiments: The reactor 
experiments do not depend on 5qp, and the wrong fit hierarchy coincides with the right hi- 
erarchy. For T2K, which is simulated with neutrino running only, there is some dependence 
on 5qp, but the correlation between Sep and 8 13 cannot be resolved without antineutrino 
running. The wrong hierarchy contours are slightly shifted, but the minimum % 2 is close to 
zero. NOz/A, on the other hand, has both neutrino and antineutrino running in our simu- 
lation, which means that the correlation can, at least in principle, be resolved. The wrong 
hierarchy can in some cases be excluded because of matter effects. In the combination of 
the experiments, the combination between Daya Bay and the beams allows for a substantial 

3 Here the "right" fit hierarchy is the same as the simulated hierarchy, and the "wrong" fit hierarchy is 
the other fit hierarchy. 
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Double Chooz T2K NOvA Daya Bay Combined 




Figure 2: Fits in the 0i3-S C p plane for sin 2 26*i 3 = 0.1 and S C p — tt/2 (upper row) and Sep = 3ir/2 
(lower row). A normal simulated hierarchy is assumed. The contours refer to la, 2a, and 3a (2 d.o.f.). 
The fit contours for the right fit hierarchy are shaded (colored), the ones for the wrong fit hierarchy fit are 
shown as curves. The best-fit values are marked by diamonds and boxes for the right and wrong hierarchy 
respectively, where the minimum y 2 for the wrong hierarchy is explicitely shown. 
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Figure 3: Same as Fig. [2 for the inverted simulated hierarchy. 
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Figure 4: Fits in sin 2 ^-Am 2 ^ plane for two different sets of simulated values: sin 2 6*23 = 0.5, sin 2 26>i3 = 
(upper row) and sin 2 6*23 = 0.4, sin 2 2#i3 =0.1, <5cp = (lower row). The upper right panel does not change 
significantly if the reactor experiments are added. In the lower right panel, the unshaded contours are for the 
combination of all experiments without reactor experiments. The currently allowed MINOS+ atmospheric 
region is shown as dashed curves in all panels (3cr confidence level) [35]. The contours corresponds to the 
la, 2cr, and 3cr confidence level, respectively (2 d.o.f.). 

reduction of the allowed parameter space due to almost orthogonal measurements. In the 
most optimistic cases, the mass hierarchy can be determined at 3<r confidence level, and 
maximal CP violation can be established at relatively modest confidence as well. However, 
note that these optimistic cases represent only a very small fraction of the parameter space. 

3.3 Leading atmospheric parameters 

No matter if 13 will be discovered or not, the upcoming beam experiments will allow for 
precision measurements of the leading atmospheric parameters by exploring the dis- 
appearance channel. We illustrate the improvement of the currently allowed parameter 
range (dashed curves) in the sin 2 ^-Amjj plane in Fig. H] for the two different scenarios 
sin 2 #23 = 0.5, sin 2 2# 13 = (upper row) and sin 2 #23 = 0.4, sin 2 29 13 = 0.1, 5 C p = (lower 
row). In the maximal mixing case, the allowed region can be substantially reduced by the 
beam experiments, especially the Am| x interval, with hardly any dependence on the reactor 
experiments or the value of # 13 . 
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Figure 5: Fit range in #23 as a function of the true #23 for different experiments (la, 2a, 3a for 1 d.o.f.). 
The upper left and right, and the lower left panels are computed for sin 2 26>i3 = 0, whereas the lower right 
panel is for sin 2 2^13 = 0.1, <5cp = 0. The qualitative picture is hardly affected by the reactor experiments 
or large #13, unless both conditions apply (large #13 and reactor experiments). 

In the lower row of Fig. HI a relatively large deviation from maximal atmospheric mixing 
is chosen. Again, the allowed parameter ranges can be significantly reduced by the beams, 
apart from an octant ambiguity in 823 [37]. As illustrated in the lower right panel this am- 
biguity might be resolved if #13 is large and beams are combined with an accurate reactor 
experiment, compare the shaded (including reactors) and unshaded (excluding reactors) re- 
gions. In this case the two solutions corresponding to the two ^-octants for the appearance 
channel of the beams are located at rather different values of #13. Therefore, an independent 
determination of t9 13 from reactors can in principle resolve the ambiguity [18,22]. 

In order to discuss the sensitivity of deviations from maximal atmospheric mixing and 
the octant resolution potential as a function of the true 823, we show in Fig. [5] the fit range 
in 823 as a function of the true 823 for different experiments. Deviations from maximal 
atmospheric mixing are theoretically interesting because they may point towards deviations 
from a symmetry, see Ref. [38] for a more detailed discussion. In Fig. [5j such a deviation 
can be established if sin 2 823 = 0.5 can be excluded at the vertical axis. Obviously, this 
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measurement does not require the presence of the reactor experiments or large 613. The 
combination of the beams can establish deviations from maximal atmospheric mixing for 
I sin 2 # 23 - 0.5| > 0.07 (3o-). 

The sensitivity to the octant of 623 can be interesting to distinguish certain models. For 
instance, sum rules such as 9 2 3 ~ 7r/4 ± 6> 2 3 /2 might be tested, where the sign depends 
on the model assumptions [39]. If deviations from maximality are due to renormalization 
group effects, in a large class of models the sign of the deviation is different for normal and 
inverted hierarchy, and depends on the presence of Supersymmetry [40]. As illustrated by 
the lower two panels in Fig. [5j both the combination with the reactor data and a large #13 
are necessary to perform this measurement with the given experiments. Sensitivity to the 
octant of 623 is present if any 823 in the wrong octant can be excluded. There is no octant 
sensitivity in the upper and lower left panels of Fig. [5j From the lower right panel, we can 
read off octant sensitivity if sin 2 #23 < 0.39 or sin 2 623 > 0.61 (3<r, sin 2 2#i3 = 0.1). 

4 Sensitivity time evolution 

In the following, we discuss the sensitivity evolution within the coming years for different 
experiments, based as much as possible on official statements of the collaborations. Although 
the assumed schedules and proton beam plans may turn out to be not realistic in some cases, 
our toy scenario will be illustrative to show the key issues for the individual experiments 
within the global neutrino oscillation context. We will show the sensitivities as a function 
of time assuming that data are continously analyzed and results are available immediately. 
The key assumptions for our toy scenario are as follows. 

Double Chooz starts 09/2009 and runs 1.5 years with far detector only, then with far and 
near detector [41]. We assume that the experiment ends after five years. 

Daya Bay Starts 07/2011. At this time, it is assumed that near and far detector halls 
are completed, and all modules are ready [42]. Furthermore, all three pairs of reactor 
cores are assumed to be online. Again, we limit the operation time to five years. 

RENO starts 06/2010 [43] with both detectors and runs for 5 years. 

T2K starts 09/2009 with virtually MW beam power, which increases linear to 0.75 MW 
reached in 12/2012. From then we assume the full target power of 0.75 MW. (This is 
an approximation for the T2K proton plan from Ref. [44]). The Super-Kamiokande 
detector is online from the very beginning and the beam runs with neutrinos only, (at 
least) until 2018 or 2019. 

NOz/A starts 08/2012 with full beam (0.7 MW), but 2.5 kt detector mass only. Then the 
detector mass increases linearly to 15 kt in 01/2014. From then we assume the full 
detector mass of 15 kt [45]. The beam runs with neutrinos first, until the equivalent 
of three years operation at nominal luminosity (cf., Tabled]) is reached, i.e., 03/2016. 
Then it switches (possibly) to antineutrinos and runs at least until 2019. 



10 



sin 2 26> 13 sensitivity limit (NH, 90% CL) 



10" 



o 

03 
> 



10" 



10 u 



GLOBES 2009 




/ CHOOZ+ 
i 

Solar excluded 



NOvA:vonly 



2010 



2012 



2014 
Year 



2016 



2018 



Figure 6: Evolution of the #i 3 sensitivity limit as a function of time (90% CL), i.e., the 90% CL limit 
which will be obtained if the true #13 is zero. 



4.1 Finding versus constraining #13 

We now discuss the sensitivity of the different experiments to #13. We include two quali- 
tatively different aspects in the discussion: The #13 sensitivity limit and the # 13 discovery 
potential. The #13 sensitivity limit describes the ability of an experiment to constrain # 13 if 
no signal is seen. It is basically determined by the worst case parameter combination which 
may fake the simulated 813 = 0. The sensitivity limit does not depend on the simulated 
hierarchy and 5qp, as the simulated 613 = 0. For a more detailed discussion, see Ref. [21], 
App. C. The 813 discovery potential is given by the smallest true value of 8 13 > which 
cannot be fitted with 9 13 = at a given CL. Since the simulated # 13 , <5 CP , and hierarchy 
determine the simulated rates, the 613 discovery potential will depend on the values of all 
these parameters chosen by nature. On the other hand, correlations and degeneracies are of 
minor importance because for the fit #13 = is used. The smallest #13 discovery potential 
for all values of 5cp and the MH (risk-minimized #13 discovery potential) is often similar to 
the 9\3 sensitivity limit. This holds to very good approximation for reactor experiments, 
where statistics are Gaussian and the oscillation physics is simple. For beam experiments 
differences occur due to Poisson statistics as well as more complicated oscillation physics 
implying correlations and degeneracies. 

We show the 613 sensitivity limit as a function of time in Fig. We observe that the 
global sensitivity limit will be dominated by reactor experiments. As soon as operational, 
Daya Bay will dominate the global limit. For Daya Bay, time is not critical, but matching 
the systematics or statistics goals is. If the assumed schedules of both, Double Chooz and 
Daya Bay are matched, Double Chooz will dominate the 613 sensitivity for about two years 
in the absence of RENO. If available, RENO, on the other hand, will dominate the 613 
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Figure 7: Evolution of the #13 discovery potential as a function of time (90% CL), i.e., the smallest value 
of 0i3 which can be distinguished from zero at 90% CL. We assume the normal and inverted simulated 
hierarchies in the left and right panels, respectively. The bands reflect the (unknown) true value of <5cp- 

sensitivity if it is operational significantly before the end of 2011. Since we do not obtain 
any CP violation or mass hierarchy sensitivity before 2014, as we shall demonstrate later, the 
reactor contribution to those will be completely dominated by Daya Bay. As a peculiarity, 
the #13 sensitivity of NOz/A is improved by switching to antineutrinos. However, the global 
limit will at that time be dominated by the reactor experiments. 

The #13 discovery potential is shown in Fig. [7J as a function of time. For the beam 
experiments, the dependence on the true value of 5cp is shown as shaded region, whereas the 
reactor experiments are not affected by the true <5cp- There is a small dependence on the true 
mass hierarchy for the beam experiments, compare left and right panels. The comparison 
of Figs. [7J and [6] shows that suitable values of 5qp may significantly improve the discovery 
potential of beams compared to their sensitivity limit. Indeed, the beam experiments may 
discover 613 for smaller #13 than Daya Bay in a small fraction of the parameter space (see also 
Fig. [1]). Overall, it may however be more likely that the reactor experiments are faster. The 
NOi^A bands become more narrow for some additional antineutrino running, which means 
that the best case potential gets slightly worse, but the worst case becomes somewhat better. 
Again, this may not be an argument for antineutrino running since the anticipated Daya 
Bay sensitivity may not be exceed-able. For a more detailed discussion of the potential 
antineutrino running, we refer to Sec. [5j 

Note that this discussion is based on the unitarity standard three-flavor oscillation frame- 
work. If the search for new physics is taken into account, different reactor experiments, or 
reactor experiments and superbeams, may imply different information and therefore be very 
complementary; see, e.g., Refs. [46,47]. 
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Figure 8: Mass hierarchy (left panels) and CP violation (right panels) discovery potentials as a function 
of true sin 2 2#i3 and fraction of true <5cp at 90% CL from T2K, NO^A and reactor data. The upper row 
corresponds to the normal simulated hierarchy, the lower row to the inverted simulated hierarchy. The 
different shadings correspond to different points of time, as marked in the plots (note that "2015" here 
means mid 2015). The dashed curves refer to NO^A with neutrino running only, whereas the shaded 
contours refer to the nominal NO^A neutrino-antineutrino plan. If no contour is shown for a specific year, 
there is no sensitivity. Note the different scales on the vertical axes. 
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4.2 Mass hierarchy and CP violation discovery potentials 

We show in Fig. [8] the mass hierarchy (left panels) and CP violation (right panels) discovery 
potentials as a function of true sin 2 29^ and fraction of (true) 5qp from T2K, NOz/A, and 
Daya Bay. The upper row corresponds to the normal simulated hierarchy, the lower row 
to the inverted simulated hierarchy. The different contours represent different points in 
time, and can be viewed as timeslices. Obviously, there will be no mass hierarchy discovery 
before 2014, and no CP violation discovery before 2015. In the most optimistic case without 
upgrades, the mass hierarchy and CP violation can be discovered in about 50% of all cases 
of #cp by 2019. Note, however, that we show the 90% CL potentials here, whereas there is 
only extremely poor 3a sensitivity. Therefore, as a first conclusion, there will be no high 
significance determination of the mass hierarchy or CP violation without the next generation 
of experiments. In the most optimistic case, some hints may be obtained in about 2015-2018. 
As discussed earlier (c/., Fig. [TJ, NOz/A plays a key role in these discovery potentials. 

As far as the NOz/A switching to antineutrinos is concerned, we show in Fig. [8] the case 
if NOz/A runs with neutrinos only as dashed curves. As expectable, for CP violation the 
antineutrino run is mandatory. For the mass hierarchy discovery, however, it is of secondary 
importance. For a more detailed discussion, see Sec. [5j 

Interesting observations can be obtained from the comparison between Fig. [8] and the 
9 13 time evolution. If sin 2 2$i 3 < 0.02, it may already be excluded by Daya Bay early 2012 
(c/., Fig. EJ). In this case, it will be clear already before the NOz/A startup that there is 
no chance to find the mass hierarchy or CP violation without significant upgrades or new 
experiments. On the other hand, if sin 2 26 l 1 3 > 0.02, this will be known at about the same 
time (cf., Fig. [7j), meaning that NOz/A has a realistic chance to see something. In this way, 
the "branching point" 2012 will be an interesting point of time at which strategic decisions 
on the future neutrino oscillation program can be made, such as in favor of upgrades or a 
new high intensity facility. 

5 Beam upgrades for large #13 and the v-i> optimization 

In the previous sections we have seen that the sensitivity to CPV and MH of the discussed 
experiments in their nominal configuration as defined in Sec. [2] is marginal, at best. In this 
section we address the following question. Imagine that a finite value of 6*13 is discovered 
soon; can "modest upgrades" of the experiments considerably improve the sensitivity to 
CPV and MH? With "modest upgrades" we mean modifications of existing equipment and 
infrastructure. This includes a longer running time and an upgraded beam power for both 
accelerator experiments and the addition of antineutrino running in T2K0 It does not 
include new beam lines or new detectors. In particular, our toy scenario is the following]! 

T2K We assume that a proton driver is installed, which increases the beam power from 
0.75 to 1.66 MW, linearly from 2015 to 2016 [49]. 

4 The anti-neutrino running beam fluxes for the 2.5 degree off-axis beam for T2K are taken from [48]. 
5 In this section, if not stated otherwise, years denote always the middle of year, i.e., 2019 = July 1st, 
2019. 
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Figure 9: The contours show the CP fraction for which CP violation can be established at 90% CL as a 
function of the year and the fraction of antineutrino running fp up to that year. We assume the upgraded 
run plan. The black contours assume T2K (right hand panel) or NO^A (left hand panel) data only, whereas 
the colored contours show the result for the case where reactor data on #13 is included. 

NOz/A At Fermilab, the proton driver "ProjectX" is discussed. We assume a linear increase 
from 0.7 to 2.3 MW from March 2018 to March 2019 [50]. 

For Double Chooz, in principle, there exists the upgrade option to Triple Chooz [29]. How- 
ever, from Fig. [6j we do not see how that could compete with Daya Bay (either timescale- 
wise, or physics- wise) . Therefore, we do not discuss this possibility here. Likewise, neither 
for RENO nor for Daya Bay upgrade options are discussed/anticipated. 

To be specific, we assume in the following that the true value of #13 is sin 2 2#i3 = 0.1 
and the true hierarchy is normal. According to Fig. [7] this value of # 13 will be discovered 
before 2011 by Double Chooz, T2K and RENO. We focus first on CPV and comment on 
the MH in Sec. 

5.1 Optimization of neutrino/antineutrino running 

We start our analysis by discussing the optimization of the fraction of neutrino/antineutrino 
exposures. Fig. [9] shows the fraction of <5cp-values for which CPV can be established at 
90% CL as a function of time and as a function of the fraction of antineutrino running 
f p for NOi^A (left) and T2K (right). Here, fp = 0(1) corresponds to the full exposure 
with neutrinos (antineutrinos). As expected, in order to address CPV a sizable fraction of 
antineutrino running is required: fp = 0.6 — 0.7 for T2K and fp = 0.5 — 0.6 for NOz/A. 
The plot is for normal hierarchy, but results are similar for inverted hierarchy. Without 
reactor data fp = is strongly disfavored even when the beams are upgraded, indicating 
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Figure 10: The contours show the CP fraction for which CP violation can be established at 90% CL as 
a function of the antineutrino fraction in NOz/A and T2K in the year 2019 assuming the nominal run plan 
(left hand panel) or in the year 2025 with the upgraded run plan (right hand panel). The black contours 
assume the combined data from NOz/A and T2K only, whereas for the colored contours also reactor data 
on 9i3 is included. The thick contours have a spacing of 0.1 whereas the thin contours have a spacing of 
0.02. The black disk denotes the optimal point for beam data only and whereas the colored disk denotes the 
optimal point for beam and reactor data combined. The circles denote the optima found for each experiment 
individually. 

that spectral information is not enough to provide a measurement of Sep- 

NOz/A with upgrades reaches a CP fraction of about 60%, while T2K achieves only 
about 30%. The reason is that the beam upgrade is a factor of 3.2 for NOzvA whereas it 
is only 2.2 for T2K. Moreover, antineutrino running is very difficult for T2K due to the 
low beam energy Since the ratio of antineutrino/neutrino cross sections as well as fluxes 
decreases with energy, statistics is rather poor for T2K antineutrino data. This is especially 
true for the 2.5° off-axis T2K beam. Optimization of beam optics or energy is beyond the 
scope of this work. 

The blue (light gray) curves in Fig.[9]show the case when reactor data on #13 are combined 
with that of the beam experiment. While for NOz/A, reactor data have only a minor impact, 
for T2K, the combination with a reactor experiment can considerably relax the requirement 
on f 9 and even replace the antineutrino running altogether in most cases [51]. However, the 
optimal CPV fraction for T2K is not improved significantly by adding reactor data. 

Fig. [TD] shows the neutrino/antineutrino optimization for the combined data of T2K 
and NOzvA, with and without reactor data. In the left panel we assume only the nominal 
luminosities for the beams, but allow for an arbitrary fraction of antineutrinos for both 
experiments. Without reactor data, the optimum is at fp° uA / fJ 2K = 0.5/0.7. The CP 
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fraction is about 32% and antineutrino data is required for both experiments. If reactor 
data are added, the optimum is at fp° uA / fp 2K = 0.5/0, which is remarkably close to the 
nominal plan for the experiments. The CP fraction is 40%, in agreement with Fig. [HI We 
observe again that reactor data are more useful than the statistically weak antineutrino 
data in T2K, which means that it is preferable to fully explore neutrino data in T2K. 

The right panel of Fig. [TO] shows the same analysis for the upgraded beams. The ten- 
dency is to use large antineutrino fractions in NOz/A and less in T2K, the optimum being 
fv° vA lfv 2K = 0.9/0.5(1/0.1) without (with) reactors. We find an optimal CP sensitivity 
at 90% CL for beam+reactor data for about 76% of all <5cp values. The landscape be- 
comes rather flat and similar CP fractions can be obtained for a wide range of antineutrino 
fractions. 



5.2 Optimal run plan 

In the previous subsection we investigated the optimal total fraction of antineutrino running 
for the full lifetime of the experiment. In this section, we proceed further and ask the 
following question: what is the neutrino/antineutrino running strategy such that at at each 
point in time the optimal sensitivity to CPV is obtained? Let us first define a reference 
neutrino-antineutrino scenario ( "nominal run plan" ) as follows: 

Phase I: Until design luminosity reached T2K: neutrinos only, NOz/A: 50% neutrinos 
followed by 50% antineutrinos (such as in LOIs) 

Phase II: After design luminosity reached until mid 2025 T2K: 50% neutrinos fol- 
lowed by 50% antineutrinos. NOz/A: 50% neutrinos followed by 50% antineutrinos 

Note that the definition of phase I and II is not connected with the upgrade point of time, 
but the design luminosity. The #cp fractions for CPV at 90% CL for this nominal run plan 
are shown in Fig. [11] as a function of time with blue bands, where light and dark colors 
indicate neutrino and antineutrino running, respectively. 

Now we investigate whether a globally optimized neutrino-antineutrino run plan has 
the potential to improve the performance with respect to this reference plan. This (non- 
deterministic) optimization is performed with the help of a genetic algorithm described in 
detail in appendix [A] Roughly, the algorithm evolves a set of switching times between 
neutrinos and antineutrinos for T2K and NOz/A by favoring "individuals" with the largest 
CP fraction at each point in time and disfavoring the ones with too many switching points. 
After evolving a randomly chosen population of 2000 individuals for 50 generations we 
average over the 100 "fittest" individuals. The results of this optimization search are shown 
as green bands in Fig. [TTJ see also Tab. [3] in the appendix. 

Let us first discuss the nominal exposures without the beam upgrades (upper row of 
panels). We find that the reference z/-z/ run plan outlined above provides a very suboptimal 
performance and can delay physics by four to five years for beam data without reactors. 
If reactor data are added, delays will be reduced to one or two years. Again we observe 
that reactor data allow to avoid antineutrino running in T2K, and we find that an early 
antineutrino run in NOz/A is preferred. The optimal strategy is not costly in terms of 
switches, only one or two additional switches are required. 
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Figure 11: Fraction of <5cp values for which CP violation can be discovered at 90% CL as a function 
of time for the nominal run plan (upper row) and the upgraded run plan (lower row). These plots assume 
the true sign of Am| : > 0. The left hand column shows the results for a combination of T2K and NO^A 
alone, whereas the right hand column includes reactor data. The green band is for the globally optimal 
run plan, whereas the blue band is for the nominal run plan. Light and dark colors indicate neutrino and 
antineutrino running, respectively. The black circles denote the absolute maximal performance for a given 
time. 
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The lower panels in Fig. [IT] show the results for upgraded beams. Again the reference 
v-v run plan provides a suboptimal performance and can delay physics by two to three 
years for beam data only. Inclusion of reactor data helps a lot and delays are reduced to 
one year. Here reactor data no longer allow for circumventing antineutrino running in T2K. 
We observe that the final overall performance with and without reactors is very similar if 
the run plan is optimized. Thus, it seems that one should strive to optimize the beams as 
much as possible on their own. In that way the reactor data can be used as an independent 
cross check of the result, for instance, to provide some sensitivity to non-standard neutrino 
physics. Again we find that the optimal strategy prefers early antineutrino data from NOz/A, 
and that only one or two additional switches are required. 

5.3 Time evolution of physics potential 

Having optimized the CPV performance for a true value sin 2 2# 13 = 0.1 in the previous 
subsections, we now relax this assumption^! and discuss the sensitivities of these optimized 
configurations as a function of #13. We consider the upgraded beams for T2K and NOz/A 
combined with reactor data, and show results for CPV as well as for the neutrino mass 
hierarchy. Fig. [12] shows the discovery potential as a function of true sin 2 2# 13 and fraction 
of true #cp for times from 2015 to 2025. The upper row of this figure shows the discovery 
potentials at the 90% CL. These results can be compared to Fig. [S] (upper panels), where one 
should note the different scales on the vertical axes. Obviously, with the optimal upgrade 
plan, there is a significant improvement of the MH and CPV discovery potentials. At the 
90% confidence level, there will be hints for the MH and CPV for sin 2 2#i3 > 0.05 for most 
values of #cp around 2025. 

However, certainly a 90% CL is not sufficient to make any meaningful statement about 
a discovery. Therefore, we show in the lower row of Fig. [12] the corresponding results at 
3a CL. Obviously the sensitivity regions reduce drastically, however, we see from the figure 
that assuming both beams upgraded, a fully optimized neutrino/antineutrino run plan, and 
data from reactors a non-negligible discovery potential at 3a will be reached in 2025. The 
mass hierarchy can be identified for sin 2 29i 3 > 0.05 for about 20% to 40% of 5qp values, 
whereas CPV can be discovered for sin 2 29 13 ^ 0.02 for 25% of 5qp values. In both cases, 
MH and CPV, there is sensitivity for values of 5qp around 37r/2 (tt/2) if the true hierarchy 
is normal (inverted). This behaviour is visible in Figs. [2] and [3] and it is related to the sign 
of the matter effect, see, e.g., Ref. [52] for a discussion. 

Let us mention that in the previous two subsections we have optimized the sensitivity to 
CPV. However, it turns out that this is also very close to the optimal performance for the MH 
determination. Furthermore, it is a feature of our optimization plan that the performance in 
the case of a (true) inverted hierarchy is similar to one for the normal hierarchy. Therefore 
it is not shown explicitely here. 

6 Note, that in general the x 2_ f unc tions will have lower values and hence be less steep for smaller #13. 
Therefore, optimization is less crucial and a wider range of running fractions is acceptable. Thus a solution 
which was optimal for large #13 will still be very close to optimal for somewhat smaller 6*13. However, all 
sensitivities will be decreased. 
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Figure 12: Mass hierarchy (left panels) and CP violation (right panels) discovery potentials as a func- 
tion of true sin 2 2# 13 and fraction of (true) <5cp for our optimal run plan including upgrades and reactor 
experiments; c/., Table [31 The upper panels are for 90% CL, the lower panels for 3<r CL. The different 
shadings corresponds to different points of time, as marked in the plots. Note the different scales on the 
vertical axes compared to Fig. [8] 
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sin 2 26»i3 |Am^| | sin 2 6> 23 - 0.5| 



Double Chooz 0.033 (0.060) 

T2K 0.004-0.027 (0.011-0.040) jjjjg 0.055 (0.074) 

RENO 0.018 (0.033) 

Daya Bay 0.007 (0.012) 

NOz/A 0.005-0.014 (0.011-0.025) + _lfj C^B) - 065 (0-092) 

Table 2: Summary of nominal sensitivities at 90% (3a). We show the discovery potential for sin 2 2#i3 
(where for the beams we give the best and worst sensitivity depending on Sep and mass hierarchy), the 
accuracy on jAm^l (for maximal mixing), and the sensitivity to deviations from maximal 623 mixing. 



6 Summary and conclusions 



In this work we have discussed the physics potential of the upcoming neutrino oscillation 
experiments Double Chooz, Daya Bay, RENO, T2K, and NOzM.. In the first part we have 
reconsidered the sensitivities from their nominal exposures, such as stated in the experimen- 
tal proposals. These results are summarized in Table El While Double Chooz and maybe 
RENO could allow for a fast discovery of a relatively large value of #13, the ultimate reactor 
experiment will be Daya Bay, which will dominate a few months after it comes online. The 
#13 performance of the beam experiments will depend on the value of 5qp and the perfor- 
mance indicator. In the case of no signal, the limit on #13 will be similar to the one from 
Double Chooz, whereas a favorable value of 5qp will allow for a discovery for slightly smaller 
#13 than Daya Bay. 

We have found that the global sensitivity to CP violation and the neutrino mass hierarchy 
of all experiments with nominal exposures is marginal. For the largest allowed values of #13, 
typically a hint at 90% CL can be obtained for about 25% to 50% of all possible values of 
Sep, while almost nothing can be said at 3a. Therefore, we have investigated the possibilities 
to increase the sensitivity by minor upgrades of the beam experiments in case #13 is not 
too far from its present bound and hence discovered soon. These upgrades are based on 
existing equipment and include an increase of beam power with the help of proton drivers, 
longer running times, and the addition of antineutrino running in T2K. We have performed 
an optimization study concerning the distribution of neutrino and antineutrino data runs 
in T2K and NOz/A in order to maximize the global sensitivity reach. We have found 
that typically communication between the two beam experiments will improve the overall 
sensitivity. For example, the optimal sensitivity usually requires relatively early antineutrino 
data from NOz/A. Furthermore, we have found that the sensitivity of the optimized running 
strategy for T2K plus NOz/A is rather similar to the one with reactor experiment data used 
in addition. 

These results are summarized in Fig. [TSJ Assuming both beams upgraded, a fully opti- 
mized neutrino/antineutrino run plan, and data from reactors, we have found that a hint 
for CP violation at 90% CL can be obtained around 2025 for a large fraction of 6qp for a 
reasonably large #13. However, the discovery potential for CP violation and mass hierarchy 
drastically reduces at the 3cr CL. Nevertheless, a non-negligible discovery potential at 3<r 
will be reached in 2025: The mass hierarchy can be identified for sin 2 2# 13 > 0.05 for about 
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20% to 40% of Sep values, whereas CP violation can be discovered for sin 2 2^13 > 0.02 
for 25% of Sep values. Sensitivity to the mass hierarchy and CP violation is obtained for 
values of S C p close to maximal CP violation: 5 C p — 3n/2 (ir/2) for a true normal (inverted) 
hierarchy, but not for the opposite case. 

Let us mention that NOz/A dominates the MH determination and typically also has a 
somewhat better sensitivity than T2K for CPV. For example, at sin 2 2# 13 = 0.1, NOz/A 
alone obtains a CP fraction for CPV at 90% CL of about 60%, whereas T2K reaches 35%, 
and the combination yields 76% (both beams with upgrades). The reasons are that the 
beam upgrade we consider is a factor of 3.2 for NOz/A whereas it is only 2.2 for T2K, and 
antineutrino running is very difficult for T2K due to the low beam energy. Since the ratios 
of antineutrino/neutrino cross sections as well as fluxes decrease with energy, statistics is 
rather poor for T2K antineutrino data. 

Our results raise the question on how to adapt the global oscillation strategy in case #13 
were discovered soon. Although "minor upgrades" of existing facilities may provide a non- 
negligible sensitivity to CP violation and the mass hierarchy, there is high risk associated 
with this strategy, since for about 75% of all possible values for 5 C p no discovery will be 
possible at the 3a level. In contrast, a high precision facility such as a wide band superbeam 
or a beta beam with very large detectors or a neutrino factory [25,26] would certainly be 
able to perform a solid determination of CP violation and the mass hierarchy in case of 
such a large value of #13. Therefore, we conclude that the upcoming generation of oscillation 
experiments may lead to interesting indications for the mass hierarchy and CP violation, but 
it is very likely that an experiment beyond the upcoming superbeams (including reasonable 
upgrades) will be required to confirm these hints. 
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A Optimization method 



Luminosity C is given by the time integrated product of target mass m and beam power p 

C(t) := f dt'm(t')p(t'), (1) 

Jto 

where we allow for m and p to be time dependent. We define the normalized luminosity L 

by 

L(t) = ^ , (2) 

where Cq corresponds to the design luminosity of the experiment. Since, in the following 
we will discuss T2K and NOz/A every quantity pertaining to T2K will carry T as subscript 
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and those carrying subscript N pertain to NOi^A. 

T2K and NOi^A both can run neutrinos or antineutrinos. We will denote the the fraction 
of antineutrinos in each experiment by /, where 



The current run plans foresee no antineutrino running for T2K, fj, = and 50% of antineu- 
trino running for NO^A, /jy = 0.5. For a given, true value of sin 2 2^13 we can compute the 
obtainable fraction of 5qp, C, for the measurement of CP violation as a functions of both 
/t and /at. This is shown in figure [101 Of course this will be a time dependent quantity de- 
noted by C(f]y(t), frit))- To simplify our notation we will write C t (fN, fr)- For algorithmic 
purposes, it turns out to be useful to define the normalized CP fraction C by 

Ct(f N jT) = CtifN /i JT) with e t = maxC t (f N ,f T ). (4) 
Ct jnJt 

Thus, Ct will have values from to 1 for all times t. Note, that there will be seperate and 
distinct C t and C t for the true hierarchy normal (NH) or inverted (IH). 

Since an experiment needs to reverse the polarity of the electric current in the focusing 
horn, its operation will be divided into phases of neutrino running and those of antineutrino 
running. Any such division can be described a set of times at which the operation mode 
switches, plus the initial running configuration. We will denote such a set as g = {gi, . . . ,gi g } 
for initial neutrino running and with g = {gi, . . . ,gi s } for initial antineutrino running. The 
gi denote the individual switching times and l g the total number of switches, g uniquely 
determines the the antineutrino fraction as a function of time 



with g k+1 <t, (5) 
= 1. Thus, we now can fix the 

K{t,g N ,g T ) = C t (f 9N (t)J gT (t)) , (6) 

in case both experiments started with neutrino running. Again, there are two functions 
k nh (t, gw, <7t) f° r true hierarchy normal and k ih (t, g^, g?) for true hierachy inverted. The 
goal now is to find a sequence of switching times g which yields n NH (t, g^, g?) = 1, Vt and 
K IH (t, <7jv, grp) = 1, Vi. In this case, we would have optimal sensitivity at any given moment 
in time, irrespective of the true mass hierarchy. We will search such g or g by using a genetic 
algorithm. To this end, we define a fitness function to evaluate the relative performance of 
different sequences g. 



_ L(t) - L(g k+1 ) ^ L(g 2i +i) - L(g 2i ) 
h[) ~ L{t) L[t) 

for starting with antineutrinos otherwise the sum starts at i 
value of C t for any t by inserting the f„{t) like this 



^(gN,gr) jqxx^i _ K NH( Th g Nj g T ^i _ K IH (Ti,g N ,g T )} l 9N + l gT ' 
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where the r< are suitably chosen points in time, e.g. every two years from to on. The 
first terms selects those solutions which manage to balance the deviation from optimal 
performance at all times. The second term penalizes solution which switch very frequently 
without gaining significantly in the first term. The genotype of each individual is given by 
a pair of {gN.gr}- In our implementation of the algorithm the length of each sequence l 9N 
and lg T are free parameters and subject to evolution. We start with a population of widely 
varying sequence length. Thus, we do not have to specify the number of switches. Also, two 
times Qi and g^\ are merged if they are less than a month apart. Each population consists 
of 2000 individuals which are evolved over 50 generations. We consider 4 populations 

{9n,9t}, {9n,9t}, {9n,9t}, {9n,9t}- (8) 

The results for each population are shown in table [31 The Tj are 

{2010.5, 2011.5, 2012.5, 2013.5, 2014.5, 2015.5, 2016.5, 2017.5, 2018.5, 2019.5} 

for nominal luminosity and 

{2010.5, 2011.5, 2012.5, 2013.5, 2014.5, 2015.5, 2017.5, 2019.5, 2021.5, 2023.5, 2025.5} 

for upgraded luminosity. 
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